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A  region  may  be  represented  by  specifying  the  curves  that  bound  it.  When  p  pro¬ 
cessors  are  available,  typical  operations  on  regions  so  represented  can  be  performed  much 
faster  than  using  a  single  processor.  This  paper  presents  parallel  algorithms  to  deter¬ 
mine  properties  of  regions,  such  as  area;  to  perform  operations  on  regions,  such  as  union 
and  intersection;  to  determine  if  a  point  lies  inside  a  region;  and  to  determine  whether  a 
given  digital  curve  could  be  the  boundary  of  a  region.  Some  of  the  algorithms  involve 
sorting,  the  time  complexity  of  which  depends  on  the  particular  model  of  parallel  compu¬ 
tation  used. 
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1.  Introduction 


A  region  in  a  digital  image  is  determined  by  specifying  the  interior  of  the 
region  or  its  boundary.  The  interior  (or  exterior)  of  a  region  can  be  represented 
by  various  schemes  such  as  run  length  codes,  medial  axis  transforms,  unions  of 
rectangles,  quadtrees,  etc.  The  boundary  of  a  region  may  be  represented  using 
chain  codes  or  crack  codes.  For  simple  regions,  these  are  usually  compact 
representations  and  allow  faster  processing.  See  [l]  for  a  detailed  discussion. 

Because  of  the  importance  of  regions  in  digital  image  processing,  and 
because  of  the  wide  availability  of  multiprocessor  architectures,  it  is  of  interest  to 
inquire  if  the  availability  of  p  processors  results  in  a  p-fold  speed-up  in  the  per¬ 
formance  of  operations  on  compact  representations  of  regions.  In  this  paper,  we 
will  consider  this  question  for  the  case  of  the  boundary  representations. 

Chain  codes  and  crack  codes  are  two  ways  to  represent  border,  in  which  a 
starting  point  and  a  sequence  of  moves  around  the  border  are  specified.  The 
chain  coding  scheme  moves  along  a  sequence  of  border  points  belonging  to  the 
region.  The  8-neighbor  chain  code  goes  from  a  point  to  one  of  its  eight  neighbors 
in  directions  that  are  multiples  of  45  * .  The  4-neighbor  chain  code  goes  from  a 
point  to  one  of  its  four  neighbors  in  directions  that  are  multiples  of  90  * .  In  the 
crack  code  representation,  the  moves  are  along  the  cracks  between  the  pixels 
belonging  to  the  region  and  the  adjacent  pixels  outside  the  region.  Each  move  is 


either  left,  right,  up  or  down.  See  Figure  1. 
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Figure  1:  Boundary  representations  of  a  region,  starting  from 
its  upper  left  corner 

Crack  code:  000303323221210121 
4-neighbor  chain  code:  00303232120112 
8-neighbor  chain  code:  007654313 
Starting  point  is  the  northwest  comer. 

Key:  In  crack  code  and  4-neighbor  chain  code,  *  =  90 i  ° ; 

in  8-neighbor  chain  code,  t  =  45i  ° .  All  angles  are  meas¬ 
ured  counterclockwise  from  the  positive  a: -axis. 


From  now  on  we  consider  only  crack  code,  and  we  use  the  letters  N,  S,  E,  W 
to  represent  the  move  segments.  Most  of  the  results  in  the  paper  also  apply  to 
4-neighbor  chain  code,  with  simple  modifications. 

To  define  the  interior  and  exterior  of  a  closed  curve,  we  use  the  implicit 
“direction”  associated  with  every  crack  code  fragment,  which  is  the  direction  in 
which  a  move  is  made.  We  adopt  the  convention  that  the  interior  of  a  closed 
curve  always  lies  to  the  right  of  the  curve  when  moving  along  the  direction 
implied  by  the  boundary  code.  With  this  convention,  we  may  also  represent 
holes  in  the  regions.  A  hole  is  represented  by  the  crack  code  of  its  boundary,  but 
since  a  hole  boundary  and  an  outer  boundary  enclose  areas  on  opposite  sides  of 
the  boundary,  a  hole  boundary  must  have  a  “direction  sense”  opposite  to  that  of 
the  outer  boundary.  We  use  the  convention  that  the  sense  of  an  outer  boundary 
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is  specified  by  a  boundary  code  in  the  clockwise  direction  while  that  of  a  hole 
boundary  is  specified  by  a  boundary  code  in  the  anticlockwise  direction.  See  Fig¬ 
ure  2. 
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hole  boundary 
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Figure  2 

The  purpose  of  this  paper  is  to  present  a  variety  of  algorithms  which  involve 
the  crack  code  representation  and  which  may  be  executed  on  suitable  multipro¬ 
cessor  networks. 

We  assume  that  we  have  only  a  fixed  number,  p ,  of  processors.  In  general 
p  «  the  number  of  boundary  points.  The  parallel  model  of  computation  we  use 
is  a  synchronous  model.  Each  of  the  p  processors  has  its  own  memory  and  its 
own  address.  At  each  time  step,  a  processor  can  send  out  a  message  specifying 
the  receiver’s  address;  selective  broadcasting  is  also  allowed.  At  each  time  step,  a 
processor  can  receive  only  one  message.  ZMOB  [2}  is  an  example  of  this  model. 
This  model  of  computation  is  strictly  weaker  than  the  shared  memory  models, 
which  are  currently  unrealizable  in  practice.  Therefore  all  the  algorithms  in  this 
paper  can  be  executed  using  any  of  the  shared  memory  models.  In  fact,  almost 
all  the  computations  in  our  algorithms  can  be  performed  efficiently  in  a  fixed 
interconnection  tree  network  with  2-way  communications  along  the  links.  The 
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only  exception  is  a  sorting  step  which  enables  us  to  find  duplicates.  A  tree  net¬ 
work  does  not  perform  this  step  very  efficiently.  Our  algorithms  can  be  executed 
efficiently  on  any  network  which  has  small  (say  log  p )  diameter  and  can  perform 
sorting  rapidly.  If  S(n  ,  p  )  denotes  the  time  required  to  sort  n  numbers  using  p 
processors,  we  can  express  the  time  required  by  our  algorithms  in  terms  of 
S(n  ,  p).  S(n  ,  p)  will  depend  on  the  model  of  parallel  computation  that  is  used. 


2.  Boundary  curves 

A  starting  point  and  a  sequence  of  (4-neighbor)  moves  can  be  regarded  as  a 
curve.  The  curve  cannot  be  the  boundary  of  a  region  unless  it  is  closed  and  does 
not  cross  itself.  In  this  section  we  present  algorithms  to  decide  if  a  given  curve 
can  be  a  region  boundary. 

First  we  describe  an  algorithm  which  obtains  the  absolute  coordinates  (with 
respect  to  the  starting  point)  of  all  points  along  the  curve.  Such  a  step  seems  to 
be  essential  in  most  of  the  computations  to  be  described,  because  the  chain  cod¬ 
ing  and  crack  coding  schemes  are  “inherently  sequential”  and  obtaining  the  abso¬ 
lute  coordinates  of  every  point  allows  parallel  processing  of  the  various  segments 
of  the  code. 

Let  n  be  the  number  of  points  in  the  code  of  some  curve.  For  simplicity  let 
2 p  +1  be  the  number  of  processors,  and  assume  p  divides  n .  We  assume  that 
the  2 p  +1  processors  can  be  connected  together  as  a  binary  tree  with  p  proces¬ 
sors  as  the  leaves.  These  p  processors  each  have  —  moves  of  the  boundary  code. 
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Algorithm  Boundary  Coordinates: 

1.  For  each  of  the  p  processors  at  the  leaves  do: 

Sequentially  compute  the  relative  coordinates  (a ,  b )  of  the  last  point  con¬ 
tained  in  that  processor  with  respect  to  the  first  point  (assumed  to  have 
coordinates  (0,0)). 

2.  For  x  :  =  0  to  log2p  -  1  do: 

Each  processor  p  at  level  i  gets  the  coordinates  (a  ,  b)  and  (c  ,  d )  stored  by 
p ’s  leftson  and  rightson  respectively,  p  stores  (a  ,  b  )  in  a  field  COUNT  and 
(a  +  c  ,  b  4-  d )  as  its  coordinates. 

3.  For  i  :  =  0  to  log2p  -  1  do: 

Each  processor  p  at  level  i  passes  the  value  it  received  from  its  father  to  its 
own  leftson.  p  adds  its  COUNT  field  to  the  value  it  received  from  its  father 
and  passes  the  sum  to  its  rightson.  The  root,  at  level  0,  passes  its  COUNT 
field  to  its  rightson  and  the  value  (0,0)  to  its  leftson. 

4.  Each  of  the  p  processors  at  the  leaves  now  sequentially  computes  the  abso¬ 
lute  coordinates  of  all  points  stored  within  it,  using  the  absolute  coordinates 
of  the  first  point  (this  is  the  value  the  processor  has  at  the  completion  of 
step  3). 

The  first  and  last  steps  each  take  time  O  |  —  j .  The  two  loops  in  steps  2,3  take 

time  O(logp)  each.  Thus  the  overall  complexity  is  +  logpj ,  where  n  is 

the  number  of  points  in  the  boundary  and  2p  +1  is  the  number  of  processors. 


It  is  clear  that  at  any  given  time  only  the  processors  at  one  level  are  active. 
Thus  it  is  possible  to  maintain  just  p  processors  at  the  leaves,  but  to  connect 
them  in  such  a  way  as  to  simulate  the  tree  connections.  Such  a  scheme,  while 
reducing  the  number  of  processors  in  half,  leads  to  a  logarithmic  increase  in  the 
maximum  degree  of  the  interconnection  network.  In  an  addressable  scheme,  such 
as  ZMOB,  or  in  a  shared  memory  model,  no  such  penalty  is  incurred.  If  we  are 
careful  to  distribute  the  boundary  codes  according  to  an  inorder  traversal  of  the 
binary  tree  with  p  processors,  using  the  same  principles  as  in  the  algorithm,  the 
above  algorithm  can  be  modified  in  the  obvious  way  without  adding  more  links 
to  the  tree  links. 

To  decide  if  a  given  curve  is  closed  or  not,  we  can  use  the  first  two  steps  of 
the  above  algorithm.  When  the  root  node  of  the  tree  obtains  the  coordinates  of 
the  last  point  in  the  curve  with  respect  to  the  first  point,  the  curve  is  closed  iff  it 

is  the  same  as  the  starting  point.  This  can  be  done  in  time  O  |  —  +  logpj . 

Next,  we  show  how  to  decide  if  a  given  chain  code  crosses  itself.  First,  the 
coordinates  of  all  the  points  on  the  curve  are  obtained  as  above.  Now,  we  need 
to  examine  these  coordinates  for  repeated  values.  One  way  to  identify  the  dupli¬ 
cate  points  is  to  sort  the  coordinates  and  examine  the  sorted  list  for  duplicate 
entries. 

In  °(f  I time  we  can  dedde  whether  this  sorted  list  contains  a  repeated 

entry.  If  this  is  not  the  case  the  curve  does  not  either  touch  or  cross  itself.  If 
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repeated  entries  are  found,  then  examine  the  directions  of  the  associated  seg¬ 
ments.  If  the  directions  of  the  two  segments  incident  on  one  point  are  E,  E  (or 
W,  W)  and  those  of  the  other  are  N,  N  (or  S,  S)  then  the  curve  crosses  itself. 
Otherwise  it  touches  itself.  Hence  the  total  time  complexity  is 

O  |  -2.  +  logp  J  +  $(n,p)  +  o|  —  |  where  S  (n  ,  p  )  is  the  time  needed  to  sort  n 

numbers  with  p  processors.  Using  the  shared  memory  model, 

5  (n  ,  p  )  =  O  [  W  fo-  n  +  2  log  p  log  n  |  [3].  We  can  determine  if  a 
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crack  code  can  be  the  boundary  of  a  region  since  a  boundary  (whether  outer  or 
hole)  must  be  closed  and  must  not  cross  itself. 

3.  Region  properties 

In  this  section,  we  will  assume  that  the  crack  code  being  considered  always 
represents  the  boundary  of  a  region  with  the  direction  conventions  assumed  in 
the  Introduction.  We  present  algorithms  to  solve  the  following  problems: 

(1)  Given  a  boundary  code,  does  it  represent  a  hole  or  an  outer  boundary? 
Determine  the  location  of  a  “corner”,  say  the  point  with  the  least  (y ,  x )  key 

value,  by  simple  comparisons  in  time  O  |—  -I-  log  pj .  Then,  observe  the 

directions  of  the  two  segments  incident  with  that  point.  If  the  direction  is 
clockwise  then  the  curve  represents  an  outer  boundary,  otherwise  it 
represents  a  hole. 


(2)  Given  a  closed  curve  and  a  point  we  can  determine  the  location  of  the  point 
relative  to  the  curve  as  follows: 

a)  Each  processor  is  given  the  coordinates  of  the  test  point.  If  broadcast¬ 
ing  is  available  this  takes  constant  time.  Otherwise,  in  a  tree  connec¬ 
tion  model,  it  takes  time  O  (log  p  ).  Next,  determine  the  coordinates  of 

all  points  along  the  curve  in  time  O  ( —  +  log  p 

b)  Each  processor  determines  in  time  0(n  /p)  if  the  test  point  coincides 
with  any  of  the  segments  of  the  boundary  contained  in  it.  If  this  is  the 
case,  the  point  is  on  the  boundary.  Otherwise  each  processor  does  the 
following: 

Determine  which  of  the  points  has  the  greatest  y  -coordinate  less  than 
the  y  -coordinate  of  the  test  point,  and  such  that: 

(i)  It  has  the  same  x  -coordinate  as  the  test  point  or 

(ii)  The  x -coordinate  of  the  test  point  falls  between  the  x- 
coordinates  of  two  adjacent  points  of  the  boundary.  A  processor 
needs  to  give  the  coordinate  of  its  last  boundary  point  to  its  next 
neighbor. 

In  either  of  these  cases,  the  y  -coordinate(s)  of  the  point(s)  are 
returned.  If  no  such  points  are  found,  a  null  value  is  returned. 

c)  Using  the  tree  connections,  determine  the  largest  y  -coordinate,  among 
all  those  passed  up,  which  is  less  than  the  y  -coordinate  of  the  test 


If  no  such  y  -coordinate  exists,  then  the  point  is  outside  the  curve. 
Otherwise,  we  observe  the  directions  of  the  segments  having  that  y- 
coordinate.  If  the  test  point  lies  to  the  right  of  the  segments,  it  is  an 
interior  point,  else  it  is  an  exterior  point.  This  step  only  takes  constant 
time. 

(3)  The  boundary  code  of  the  complement  of  the  area  enclosed  by  a  given  set  of 
boundary  codes  may  be  obtained  as  follows: 

We  assume  that  the  given  boundary  codes  are  themselves  properly  enclosed 
in  a  frame,  which  is  itself  chain  coded  as  an  outer  boundary.  The  comple¬ 
ment  can  be  obtained  by  reversing  the  direction  of  all  segments  in  the  given 
set  of  codes  in  time  O  ( n  /p  ). 


(4)  Given  a  boundary,  to  determine  the  area  enclosed  by  the  curve  we  first 
obtain  the  absolute  coordinates  of  all  points.  Then  for  each  segment  in  the 
boundary,  the  following  computation  is  performed.  If  the  segment  is  vertical 
(i.e.,  has  direction  N  or  S),  it  is  ignored.  Otherwise,  if  the  segment  has 
direction  E  (W)  and  its  y  -coordinate  is  y,  the  processor  containing  the  seg¬ 
ment  adds  (subtracts)  y  to  an  aggregate.  All  such  computations  can  be 
completed  in  time  O  ( n  Jp  ).  The  aggregates  in  the  p  processors  can  then  be 
summed  together  in  time  O  (log  p  ),  giving  an  overall  time  complexity  of 


O  |  log  p  +  —  | .  The  computation  is  equivalent  to  adding  each  of  the  areas 


shown  hatched  vertically,  and  subtracting  each  area  hatched  horizontally. 


4.  Operations  on  pairs  of  regions 


In  this  section  we  consider  tasks  involving  two  crack  codes,  both  of  which 
represent  closed  curves. 

(1)  Relative  position  of  two  regions: 

Given  two  boundary  codes  of  lengths  m  and  n  ,  and  the  absolute  coordi¬ 
nates  of  their  origins,  we  can  determine  if  the  regions  (i.e.,  the  curves)  inter¬ 
sect,  are  touching  (i.e.,  the  curves  are  tangential),  or  are  disjoint,  and  also  if 
one  region  lies  completely  within  the  other. 

We  first  obtain  the  absolute  coordinates  of  all  points  in  both  curves  in 
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time  O - H  log  p  .  The  points  appearing  on  both  boundaries  are 
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found  as  was  done  in  the  case  of  just  one  curve,  except  that  when  we  exam¬ 
ine  two  repeated  entries,  we  consider  only  points  not  both  belonging  to  the 
same  curve.  This  can  be  done  by  sorting  the  m  -f  n  points  on  the  key 
(x ,  y ).  If  the  two  curves  are  disjoint  we  can  determine  if  one  lies  inside  the 
other  by  picking  an  arbitrary  point  of  one  curve  and  solving  the  problem  of 
determining  if  this  point  lies  inside  the  other  curve,  and  if  not,  we  try  the 
same  problem  with  the  curves  exchanged. 

(2)  Given  two  boundary  codes,  to  determine  the  boundary  codes  of  the  union 
and  intersection  of  the  areas  enclosed  by  them  we  have  basically  two  cases 
to  consider:  Either  the  curves  are  disjoint,  or  they  share  points  in  common. 
This  can  be  determined  as  discussed  earlier. 

Assume  first  that  the  curves  are  disjoint.  Then,  one  may  either  lie 
inside  the  other,  or  not.  If  one  does  lie  inside  the  other: 

a)  If  both  curves  are  outer  boundaries:  the  intersection  is  the  inner 
one,  while  the  union  is  the  outer  one. 

b)  If  both  are  holes:  the  intersection  is  the  outer  one,  while  the  union 
is  the  inner  one. 

c)  If  one  is  a  hole  and  the  other  an  outer  boundary:  if  the  inner  one  is 
the  hole  boundary  nothing  need  be  done  for  the  intersection.  Other¬ 
wise  the  intersection  is  empty.  If  the  inner  one  is  the  hole  boundary, 
the  union  is  the  hole  boundary,  otherwise  nothing  need  be  done  for 
the  union. 
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We  now  consider  the  cases  when  the  curves  share  common  points. 


First,  we  determine  all  points  common  to  both  curves  by  the  procedures 
described  previously.  Each  such  point  can  find  the  address  of  the  processor(s)  it 
appears  in.  Each  processor  containing  a  shared  point  now  executes  the  following 
code: 

a)  If  the  directions  to  the  next  point  in  the  two  curves  are  opposite,  then  ignore 
the  point. 


b) 


c) 


♦ 

c 


If  the  two  curves  coincide  at  that  point  and  go  to  the  same  next  point  (i.e., 
the  next  move  directions  from  that  point  are  the  same  in  both  curves),  then 
nothing  need  be  done. 

If  the  next  moves  are  different  in  the  two  curves,  then  we  need  to  determine 
the  direction  of  the  boundary  code  from  this  point.  We  take  either  the 
intersection  or  the  union  of  the  regions  enclosed  by  the  segments  and  pick 
the  appropriate  one.  Three  basic  cases  arise: 

Case  (1):  If  A  is  a  common  point,  and  in  one  curve,  the  next  move  from 
A  is  W,  while  in  the  other  it  is  S:  the  region  enclosed  by  the 
— — «A  curve  AC  belongs  to  includes  the  area  above  AC ,  while  that  cf 

AB  includes  the  area  to  the  left  of  AB .  The  intersection  of  the 

Ib 

regions  is  the  one  included  by  AC  while  the  union  is  the  one 
included  by  AB .  Thus,  C  is  the  next  point  if  the  intersection 
was  desired,  while  B  is  the  next  point  if  the  union  was  desired. 
Once  the  successor  is  picked,  a  flag  is  put  at  the  intersection 


12 


point  to  indicate  that  the  resulting  curve  has  a  successor  in 
another  processor,  the  successor  of  this  point  is  also  stored. 

Case  (2):  Here,  the  two  curves  abut  at  A  .  In  this  case  both  B  and  C  are 
successors  of  A  on  the  two  curves.  Analyzing  the  areas  belong¬ 
ing  to  the  regions,  we  find  that  the  processor  containing  A  on 
curve  1  needs  to  make  A 's  successor  be  A ’s  successor  on  curve  2 
and  the  processor  containing  A  on  curve  2  needs  to  make  A ’s 
successor  be  A ’s  successor  on  curve  1  for  both  the  union  and 
intersection  operations.  This  is  accomplished  by  setting  flags  and 
address  as  in  case  (1). 


Case  (3):  In  this  case  the  curves  abut  along  AD .  If  the  intersection  is 


desired,  A ’s  successor  is  B  while  if  the  union  is  desired,  A ’s  sue 
cessor  is  C . 


Once  the  successors  are  picked,  a  sequential  tour  through  the  links,  starting 


from  each  intersection  point,  gives  the  boundary  code  of  the  union  or  intersec¬ 
tion. 

Once  the  common  points  are  obtained  in  time  0\ - h  log  p  +  S  (n  ,  p  )  I , 

l  p  I 


the  links  which  determine  the  union  or  intersection  can  be  found  in  time  O 


n 

P 


5.  Concluding  remarks 


WS 


*»S« 


$ 


m 


m 


Chain  codes  and  crack  codes  are  compact  representation  for  many  types  of 


regions,  in  particular,  for  simply  connected  regions.  In  this  paper,  we  have  shown 


that  using  p  processors,  each  processor  can  handle  a  piece  of  the  boundary  and 


then  combine  the  results.  This  can  be  done  quite  efficiently  as  long  as  the  proces¬ 


sors  can  be  interconnected  so  that  the  p  processors  are  not  at  more  than  distance 


log  p  apart — for  example,  as  in  a  tree  network.  In  some  of  the  algorithms,  we 


need  to  find  the  intersection  points,  and  the  sorting  speed  becomes  an  important 


factor  in  the  time  complexity. 


From  the  boundary  code  representation,  it  is  easy  to  convert  to  the  run 


length  code,  which  is  a  compact  representation  of  the  interior  of  a  curve.  Given 


the  boundary  code,  we  obtain  the  coordinates  of  all  the  points  on  the  boundary 


and  sort  them  on  the  key  (y ,  x ).  Then,  the  number  of  l’s  at  a  given  y  -level  is 


just  the  difference  between  the  greatest  and  least  x ’s  having  that  particular  y  in 


the  sorted  list.  If  the  region  contains  holes,  then  the  same  type  of  difference  gives 


the  length  of  the  hole  at  that  y  -level. 


Another  kind  of  region  representation  is  based  on  specifying  the  vertices  of  a 


boundary  polygon,  where  the  edges  are  not  constrained  to  be  rectilinear.  If  the 


intersection  points  between  two  polygonal  regions  were  to  be  determined  by  some 


other  method,  we  can  determine  the  boundary  codes  of  the  union  and  intersec¬ 


tion  by  the  same  methods  described  here. 
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